Plastomes in the holoparasitic family Balanophoraceae: Extremely high AT content, severe gene content reduction, and two independent genetic code changes.

The transition to a heterotrophic lifestyle in angiosperms is characterized by convergent evolutionary changes. Plastid genome remodeling includes dramatic functional and physical reductions with the highest degrees observed in fully heterotrophic plants. Genes related to photosynthesis are generally absent or pseudogenized, while a few genes related to other metabolic processes that take place within the plastid are almost invariably maintained. The family Balanophoraceae consists of root holoparasites that present reduced plastid genomes with an extraordinarily elevated AT content and the single genetic code change ever documented in land plant plastomes (the stop codon TAG now codes for tryptophan). Here, we studied the plastomes of Lophophytum leandri and Ombrophytum subterraneum (Balanophoraceae) that showed the remarkable absence of the gene trnE, a highly biased nucleotide composition, and an independent genetic code change (the standard stop codon TGA codes for tryptophan). This is the second genetic code change identified in land plant plastomes. Analysis of the transcriptome of Lophophytum indicated that the entire C5 pathway typical of plants is conserved despite the lack of trnE in its plastome. A hypothetical model of plastome evolution in the Balanophoraceae is presented.

Maximizing transcription of nucleic acids with efficient T7 promoters.

In vitro transcription using T7 bacteriophage polymerase is widely used in molecular biology. Here, we use 5'RACE-Seq to screen a randomized initially transcribed region of the T7 promoter for cross-talk with transcriptional activity. We reveal that sequences from position +4 to +8 downstream of the transcription start site affect T7 promoter activity over a 5-fold range, and identify promoter variants with significantly enhanced transcriptional output that increase the yield of in vitro transcription reactions across a wide range of template concentrations. We furthermore introduce CEL-Seq+ , which uses an optimized T7 promoter to amplify cDNA for single-cell RNA-Sequencing. CEL-Seq+ facilitates scRNA-Seq library preparation, and substantially increases library complexity and the number of expressed genes detected per cell, highlighting a particular value of optimized T7 promoters in bioanalytical applications.

Method to Overcome Inefficiencies in Site-Directed Mutagenesis of A/T-Rich DNA.

Site-directed mutagenesis (SDM) is an invaluable technique that enables the manipulation of DNA and therefore the primary structure and function of any encoded gene products. Commercial protocols for SDM have been optimized for Escherichia coli and mean A/T content but may hinder generation of desired products using other templates. Mutagenesis of A/T-rich DNA is often hindered by low oligodeoxynucleotide (oligo)-annealing temperatures, requiring oligos longer than manufacturer protocol recommendations. However, longer oligos can result in primer dimer formation and decreased SDM efficiencies. Commercially available kits proved inefficient at generating AT-rich mutants. We sought to generate a modified protocol that generated SDM products detectable using gel electrophoresis and that did not require an apparent limit on oligo length.

Downstream sequence-dependent RNA cleavage and pausing by RNA polymerase I.

The sequence of the DNA template has long been thought to influence the rate of transcription by DNA-dependent RNA polymerases, but the influence of DNA sequence on transcription elongation properties of eukaryotic RNA polymerase I (Pol I) from Saccharomyces cerevisiae has not been defined. In this study, we observe changes in dinucleotide production, transcription elongation complex stability, and Pol I pausing in vitro in response to downstream DNA. In vitro studies demonstrate that AT-rich downstream DNA enhances pausing by Pol I and inhibits Pol I nucleolytic cleavage activity. Analysis of Pol I native elongating transcript sequencing data in Saccharomyces cerevisiae suggests that these downstream sequence elements influence Pol I in vivo Native elongating transcript sequencing studies reveal that Pol I occupancy increases as downstream AT content increases and decreases as downstream GC content increases. Collectively, these data demonstrate that the downstream DNA sequence directly impacts the kinetics of transcription elongation prior to the sequence entering the active site of Pol I both in vivo and in vitro.

Repeat elements organise 3D genome structure and mediate transcription in the filamentous fungus Epichloë festucae.

Structural features of genomes, including the three-dimensional arrangement of DNA in the nucleus, are increasingly seen as key contributors to the regulation of gene expression. However, studies on how genome structure and nuclear organisation influence transcription have so far been limited to a handful of model species. This narrow focus limits our ability to draw general conclusions about the ways in which three-dimensional structures are encoded, and to integrate information from three-dimensional data to address a broader gamut of biological questions. Here, we generate a complete and gapless genome sequence for the filamentous fungus, Epichloë festucae. We use Hi-C data to examine the three-dimensional organisation of the genome, and RNA-seq data to investigate how Epichloë genome structure contributes to the suite of transcriptional changes needed to maintain symbiotic relationships with the grass host. Our results reveal a genome in which very repeat-rich blocks of DNA with discrete boundaries are interspersed by gene-rich sequences that are almost repeat-free. In contrast to other species reported to date, the three-dimensional structure of the genome is anchored by these repeat blocks, which act to isolate transcription in neighbouring gene-rich regions. Genes that are differentially expressed in planta are enriched near the boundaries of these repeat-rich blocks, suggesting that their three-dimensional orientation partly encodes and regulates the symbiotic relationship formed by this organism.

Ml proteins from Mesorhizobium loti and MucR from Brucella abortus: an AT-rich core DNA-target site and oligomerization ability.

Mesorhizobium loti contains ten genes coding for proteins sharing high amino acid sequence identity with members of the Ros/MucR transcription factor family. Five of these Ros/MucR family members from Mesorhizobium loti (Ml proteins) have been recently structurally and functionally characterized demonstrating that Ml proteins are DNA-binding proteins. However, the DNA-binding studies were performed using the Ros DNA-binding site with the Ml proteins. Currently, there is no evidence as to when the Ml proteins are expressed during the Mesorhizobium lo ti life cycle as well as no information concerning their natural DNA-binding site. In this study, we examine the ml genes expression profile in Mesorhizobium loti and show that ml1, ml2, ml3 and ml5 are expressed during planktonic growth and in biofilms. DNA-binding experiments show that the Ml proteins studied bind a conserved AT-rich site in the promoter region of the exoY gene from Mesorhizobium loti and that the proteins make important contacts with the minor groove of DNA. Moreover, we demonstrate that the Ml proteins studied form higher-order oligomers through their N-terminal region and that the same AT-rich site is recognized by MucR from Brucella abortus using a similar mechanism involving contacts with the minor groove of DNA and oligomerization.

HIV-1 tolerates changes in A-count in a small segment of the pol gene.

BACKGROUND: The HIV-1 RNA genome has a biased nucleotide composition with a surplus of As. Several hypotheses have been put forward to explain this striking phenomenon, but the A-count of the HIV-1 genome has thus far not been systematically manipulated. The reason for this reservation is the likelihood that known and unknown sequence motifs will be affected by such a massive mutational approach, thus resulting in replication-impaired virus mutants. We present the first attempt to increase and decrease the A-count in a relatively small polymerase (pol) gene segment of HIV-1 RNA. RESULTS: To minimize the mutational impact, a new mutational approach was developed that is inspired by natural sequence variation as present in HIV-1 isolates. This phylogeny-instructed mutagenesis allowed us to create replication-competent HIV-1 mutants with a significantly increased or decreased local A-count. The local A-count of the wild-type (wt) virus (40.2%) was further increased to 46.9% or reduced to 31.7 and 26.3%. These HIV-1 variants replicate efficiently in vitro, despite the fact that the pol changes cause a quite profound move in HIV-SIV sequence space. CONCLUSIONS: Extrapolating these results to the complete 9 kb RNA genome, we may cautiously suggest that the A-rich signature does not have to be maintained. This survey also provided clues that silent codon changes, in particular from G-to-A, determine the subtype-specific sequence signatures.

DNA helix: the importance of being AT-rich.

The AT-rich DNA is mostly associated with condensed chromatin, whereas the GC-rich sequence is preferably located in the dispersed chromatin. The AT-rich genes are prone to be tissue-specific (silenced in most tissues), while the GC-rich genes tend to be housekeeping (expressed in many tissues). This paper reports another important property of DNA base composition, which can affect repertoire of genes with high AT content. The GC-rich sequence is more liable to mutation. We found that Spearman correlation between human gene GC content and mutation probability is above 0.9. The change of base composition even in synonymous sites affects mutation probability of nonsynonymous sites and thus of encoded proteins. There is a unique type of housekeeping genes, which are especially unsafe when prone to mutation. Natural selection which usually removes deleterious mutations, in the case of these genes only increases the hazard because it can descend to suborganismal (cellular) level. These are cell cycle-related genes. In accordance with the proposed concept, they have low GC content of synonymous sites (despite them being housekeeping). The gene-centred protein interaction enrichment analysis (PIEA) showed the core clusters of genes whose interactants are modularly enriched in genes with AT-rich synonymous codons. This interconnected network is involved in double-strand break repair, DNA integrity checkpoints and chromosome pairing at mitosis. The damage of these genes results in genome and chromosome instability leading to cancer and other 'error catastrophes'. Reducing the nonsynonymous mutations, the usage of AT-rich synonymous codons can decrease probability of cancer by above 20-fold.

Horizontally acquired AT-rich genes in Escherichia coli cause toxicity by sequestering RNA polymerase.

Horizontal gene transfer permits rapid dissemination of genetic elements between individuals in bacterial populations. Transmitted DNA sequences may encode favourable traits. However, if the acquired DNA has an atypical base composition, it can reduce host fitness. Consequently, bacteria have evolved strategies to minimize the harmful effects of foreign genes. Most notably, xenogeneic silencing proteins bind incoming DNA that has a higher AT content than the host genome. An enduring question has been why such sequences are deleterious. Here, we showed that the toxicity of AT-rich DNA in Escherichia coli frequently results from constitutive transcription initiation within the coding regions of genes. Left unchecked, this causes titration of RNA polymerase and a global downshift in host gene expression. Accordingly, a mutation in RNA polymerase that diminished the impact of AT-rich DNA on host fitness reduced transcription from constitutive, but not activator-dependent, promoters.

Improved genome sequencing using an engineered transposase.

BACKGROUND: Next-generation sequencing (NGS) has transformed genomic research by reducing turnaround time and cost. However, no major breakthrough has been made in the upstream library preparation methods until the transposase-based Nextera method was invented. Nextera combines DNA fragmentation and barcoding in a single tube reaction and therefore enables a very fast workflow to sequencing-ready DNA libraries within a couple of hours. When compared to the traditional ligation-based methods, transposed-based Nextera has a slight insertion bias. RESULTS: Here we present the discovery of a mutant transposase (Tn5-059) with a lowered GC insertion bias through protein engineering. We demonstrate Tn5-059 reduces AT dropout and increases uniformity of genome coverage in both bacterial genomes and human genome. We also observe higher library diversity generated by Tn5-059 when compared to Nextera v2 for human exomes, which leads to less sequencing and lower cost per genome. In addition, when used for human exomes, Tn5-059 delivers consistent library insert size over a range of input DNA, allowing up to a tenfold variance from the 50 ng input recommendation. CONCLUSIONS: Enhanced DNA input tolerance of Tn5-059 can translate to flexibility and robustness of workflow. DNA input tolerance together with superior uniformity of coverage and lower AT dropouts extend the applications of transposase based library preps. We discuss possible mechanisms of improvements in Tn5-059, and potential advantages of using the new mutant in varieties of applications including microbiome sequencing and chromatin profiling.

The Plastid Genome of Polytoma uvella Is the Largest Known among Colorless Algae and Plants and Reflects Contrasting Evolutionary Paths to Nonphotosynthetic Lifestyles.

The loss of photosynthesis is frequently associated with parasitic or pathogenic lifestyles, but it also can occur in free-living, plastid-bearing lineages. A common consequence of becoming nonphotosynthetic is the reduction in size and gene content of the plastid genome. In exceptional circumstances, it can even result in the complete loss of the plastid DNA (ptDNA) and its associated gene expression system, as reported recently in several lineages, including the nonphotosynthetic green algal genus Polytomella Closely related to Polytomella is the polyphyletic genus Polytoma, the members of which lost photosynthesis independently of Polytomella Species from both genera are free-living organisms that contain nonphotosynthetic plastids, but unlike Polytomella, Polytoma members have retained a genome in their colorless plastid. Here, we present the plastid genome of Polytoma uvella: to our knowledge, the first report of ptDNA from a nonphotosynthetic chlamydomonadalean alga. The P. uvella ptDNA contains 25 protein-coding genes, most of which are related to gene expression and none are connected to photosynthesis. However, despite its reduced coding capacity, the P. uvella ptDNA is inflated with short repeats and is tens of kilobases larger than the ptDNAs of its closest known photosynthetic relatives, Chlamydomonas leiostraca and Chlamydomonas applanata In fact, at approximately 230 kb, the ptDNA of P. uvella represents the largest plastid genome currently reported from a nonphotosynthetic alga or plant. Overall, the P. uvella and Polytomella plastid genomes reveal two very different evolutionary paths following the loss of photosynthesis: expansion and complete deletion, respectively. We hypothesize that recombination-based DNA-repair mechanisms are at least partially responsible for the different evolutionary outcomes observed in such closely related nonphotosynthetic algae.

The Impact of Gene Silencing on Horizontal Gene Transfer and Bacterial Evolution.

The H-NS family of DNA-binding proteins is the subject of intense study due to its important roles in the regulation of horizontally acquired genes critical for virulence, antibiotic resistance, and metabolism. Xenogeneic silencing proteins, typified by the H-NS protein of Escherichia coli, specifically target and downregulate expression from AT-rich genes by selectively recognizing specific structural features unique to the AT-rich minor groove. In doing so, these proteins facilitate bacterial evolution; enabling these cells to engage in horizontal gene transfer while buffering potential any detrimental fitness consequences that may result from it. Xenogeneic silencing and counter-silencing explain how bacterial cells can evolve effective gene regulatory strategies in the face of rampant gene gain and loss and it has extended our understanding of bacterial gene regulation beyond the classic operon model. Here we review the structures and mechanisms of xenogeneic silencers as well as their impact on bacterial evolution. Several H-NS-like proteins appear to play a role in facilitating gene transfer by other mechanisms including by regulating transposition, conjugation, and participating in the activation of virulence loci like the locus of enterocyte effacement pathogenicity island of pathogenic strains of E. coli. Evidence suggests that the critical determinants that dictate whether an H-NS-like protein will be a silencer or will perform a different function do not lie in the DNA-binding domain but, rather, in the domains that control oligomerization. This suggests that H-NS-like proteins are transcription factors that both recognize and alter the shape of DNA to exert specific effects that include but are not limited to gene silencing.

Cell cycle progression in Caulobacter requires a nucleoid-associated protein with high AT sequence recognition.

Faithful cell cycle progression in the dimorphic bacterium Caulobacter crescentus requires spatiotemporal regulation of gene expression and cell pole differentiation. We discovered an essential DNA-associated protein, GapR, that is required for Caulobacter growth and asymmetric division. GapR interacts with adenine and thymine (AT)-rich chromosomal loci, associates with the promoter regions of cell cycle-regulated genes, and shares hundreds of recognition sites in common with known master regulators of cell cycle-dependent gene expression. GapR target loci are especially enriched in binding sites for the transcription factors GcrA and CtrA and overlap with nearly all of the binding sites for MucR1, a regulator that controls the establishment of swarmer cell fate. Despite constitutive synthesis, GapR accumulates preferentially in the swarmer compartment of the predivisional cell. Homologs of GapR, which are ubiquitous among the α-proteobacteria and are encoded on multiple bacteriophage genomes, also accumulate in the predivisional cell swarmer compartment when expressed in Caulobacter The Escherichia coli nucleoid-associated protein H-NS, like GapR, selectively associates with AT-rich DNA, yet it does not localize preferentially to the swarmer compartment when expressed exogenously in Caulobacter, suggesting that recognition of AT-rich DNA is not sufficient for the asymmetric accumulation of GapR. Further, GapR does not silence the expression of H-NS target genes when expressed in E. coli, suggesting that GapR and H-NS have distinct functions. We propose that Caulobacter has co-opted a nucleoid-associated protein with high AT recognition to serve as a mediator of cell cycle progression.

An AT-barrier mechanically controls DNA reannealing under tension.

Regulation of genomic activity occurs through the manipulation of DNA by competent mechanoenzymes. Force-clamp optical tweezers that allow the structural dynamics of the DNA molecule to be measured were used here to investigate the kinetics of mechanically-driven strand reannealing. When the force on the torsionally unconstrained λ-phage DNA is decreased stepwise from above to below the overstretching transition, reannealing occurs via discrete shortening steps separated by exponentially distributed time intervals. Kinetic analysis reveals a transition barrier 0.58 nm along the reaction coordinate and an average reannealing-step size of ∼750 bp, consistent with the average bp interval separating segments of more than 10 consecutive AT bases. In an AT-rich DNA construct, in which the distance between segments of more than 10 consecutive AT is reduced to ∼210 bps, the reannealing step reduces accordingly without changes in the position of the transition barrier. Thus, the transition barrier for reannealing is determined by the presence of segments of more than 10 consecutive AT bps independent of changes in sequence composition, while the length of the reannealing strand changes according to the distance between poly-AT segments at least 10 bps long.

OcculterCut: A Comprehensive Survey of AT-Rich Regions in Fungal Genomes.

We present a novel method to measure the local GC-content bias in genomes and a survey of published fungal species. The method, enacted as "OcculterCut" (https://sourceforge.net/projects/occultercut, last accessed April 30, 2016), identified species containing distinct AT-rich regions. In most fungal taxa, AT-rich regions are a signature of repeat-induced point mutation (RIP), which targets repetitive DNA and decreases GC-content though the conversion of cytosine to thymine bases. RIP has in turn been identified as a driver of fungal genome evolution, as RIP mutations can also occur in single-copy genes neighboring repeat-rich regions. Over time RIP perpetuates "two speeds" of gene evolution in the GC-equilibrated and AT-rich regions of fungal genomes. In this study, genomes showing evidence of this process are found to be common, particularly among the Pezizomycotina. Further analysis highlighted differences in amino acid composition and putative functions of genes from these regions, supporting the hypothesis that these regions play an important role in fungal evolution. OcculterCut can also be used to identify genes undergoing RIP-assisted diversifying selection, such as small, secreted effector proteins that mediate host-microbe disease interactions.

Mitochondrial genomes of praying mantises (Dictyoptera, Mantodea): rearrangement, duplication, and reassignment of tRNA genes.

Insect mitochondrial genomes (mitogenomes) contain a conserved set of 37 genes for an extensive diversity of lineages. Previously reported dictyopteran mitogenomes share this conserved mitochondrial gene arrangement, although surprisingly little is known about the mitogenome of Mantodea. We sequenced eight mantodean mitogenomes including the first representatives of two families: Hymenopodidae and Liturgusidae. Only two of these genomes retain the typical insect gene arrangement. In three Liturgusidae species, the trnM genes have translocated. Four species of mantis (Creobroter gemmata, Mantis religiosa, Statilia sp., and Theopompa sp.-HN) have multiple identical tandem duplication of trnR, and Statilia sp. additionally includes five extra duplicate trnW. These extra trnR and trnW in Statilia sp. are erratically arranged and form another novel gene order. Interestingly, the extra trnW is converted from trnR by the process of point mutation at anticodon, which is the first case of tRNA reassignment for an insect. Furthermore, no significant differences were observed amongst mantodean mitogenomes with variable copies of tRNA according to comparative analysis of codon usage. Combined with phylogenetic analysis, the characteristics of tRNA only possess limited phylogenetic information in this research. Nevertheless, these features of gene rearrangement, duplication, and reassignment provide valuable information toward understanding mitogenome evolution in insects.

Chromosome Banding in Amphibia. XXXIII. Demonstration of 5-Methylcytosine-Rich Heterochromatin in Anura.

An experimental approach using monoclonal anti-5-methylcytosine (5-MeC) antibodies and indirect immunofluorescence was elaborated for detecting 5-MeC-rich chromosome regions in anuran chromosomes. This technique was applied to mitotic metaphases of 6 neotropical frog species belonging to 6 genera and 4 families. The hypermethylation patterns were compared with a variety of banding patterns obtained by conventional banding techniques. The hypermethylated DNA sequences are species-specific and located exclusively in constitutive heterochromatin. They are found in centromeric, pericentromeric, telomeric, and interstitial positions of the chromosomes and adjacent to nucleolus organizer regions. 5-MeC-rich DNA sequences can be embedded both in AT- and GC-rich repetitive DNA. The experimental parameters that have major influence on the reproducibility and quality of the anti-5-MeC antibody labeling are discussed.

Sp7/Osterix Is Restricted to Bone-Forming Vertebrates where It Acts as a Dlx Co-factor in Osteoblast Specification.

In extant species, bone formation is restricted to vertebrate species. Sp7/Osterix is a key transcriptional determinant of bone-secreting osteoblasts. We performed Sp7 chromatin immunoprecipitation sequencing analysis identifying a large set of predicted osteoblast enhancers and validated a subset of these in cell culture and transgenic mouse assays. Sp family members bind GC-rich target sequences through their zinc finger domain. Several lines of evidence suggest that Sp7 acts differently, engaging osteoblast targets in Dlx-containing regulatory complexes bound to AT-rich motifs. Amino acid differences in the Sp7 zinc finger domain reduce Sp7's affinity for the Sp family consensus GC-box target; Dlx5 binding maps to this domain of Sp7. The data support a model in which Dlx recruitment of Sp7 to osteoblast enhancers underlies Sp7-directed osteoblast specification. Because an Sp7-like zinc finger variant is restricted to vertebrates, the emergence of an Sp7 member within the Sp family was likely closely coupled to the evolution of bone-forming vertebrates.

Single-Molecule Imaging Reveals a Collapsed Conformational State for DNA-Bound Cohesin.

Cohesin is essential for the hierarchical organization of the eukaryotic genome and plays key roles in many aspects of chromosome biology. The conformation of cohesin bound to DNA remains poorly defined, leaving crucial gaps in our understanding of how cohesin fulfills its biological functions. Here, we use single-molecule microscopy to directly observe the dynamic and functional characteristics of cohesin bound to DNA. We show that cohesin can undergo rapid one-dimensional (1D) diffusion along DNA, but individual nucleosomes, nucleosome arrays, and other protein obstacles significantly restrict its mobility. Furthermore, we demonstrate that DNA motor proteins can readily push cohesin along DNA, but they cannot pass through the interior of the cohesin ring. Together, our results reveal that DNA-bound cohesin has a central pore that is substantially smaller than anticipated. These findings have direct implications for understanding how cohesin and other SMC proteins interact with and distribute along chromatin.

Reliable handling of highly A/T-rich genomic DNA for efficient generation of knockin strains of Dictyostelium discoideum.

BACKGROUND: Social amoeba, Dictyostelium discoideum, is a well-established model organism for studying cellular physiology and developmental pattern formation. Its haploid genome facilitates functional analysis of genes by a single round of mutagenesis including targeted disruption. Although the efficient generation of knockout strains based on an intrinsically high homologous recombination rate has been demonstrated, successful reports for knockin strains have been limited. As social amoeba has an exceptionally high adenine and thymine (A/T)-content, conventional plasmid-based vector construction has been constrained due to deleterious deletion in E. coli. RESULTS: We describe here a simple and efficient strategy to construct GFP-knockin cassettes by using a linear DNA cloning vector derived from N15 bacteriophage. This allows reliable handling of DNA fragments whose A/T-content may be as high as 85 %, and which cannot be cloned into a circular plasmid. By optimizing the length of recombination arms, we successfully generate GFP-knockin strains for five genes involved in cAMP signalling, including a triple-colour knockin strain. CONCLUSIONS: This robust strategy would be useful in handling DNA fragments with biased A/T-contents such as the genome of lower organisms and the promoter/terminator regions of higher organisms.